GENERATING FUNCTIONS OF 
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Abstract. In this short note, we give simple proofs of several results and 
conjectures formulated by Stolarsky and Tran concerning generating functions 
of some families of Chebyshev-likc polynomials. 



1. Introduction 

Stolarsky [5] observed that the discriminant of the product of the polynomials 
K s (x, y) = (l + yf s + xy s and f m (y) = (y 2m+1 - l)/(y - 1) 
has a very nice factorization: 

(1.1) A y (KJ m ) = C^x 2 '-\x + 2 2s )H^(x)\ 

with Cm a rational constant and Hm \x) the following polynomial in Q[x]: 



( k-K 
\2m+ 1 



(1.2) H$(x) = Y[ ( x + 4 s cos 2s 

fc=i ^ 

This was obtained empirically for specific values of m and s. Moreover, Stolarsky 

(s) 

conjectured that the generating functions of the polynomials H m for fixed positive 
integer s are rational, with explicit empirical formulas for s = 1,2,3. In [6], Tran 
proved Eq. (1.1) with the explicit value Cm = (— l) m (2m+ l) 2m ~ 1 s 2s and obtained 
the generating function of the sequences of polynomials Hm\x) and Hm\x). 

Indeed, for 8 = 1, the polynomials Hm\x) are related to the classical Chebyshev 
polynomials of the second kind U n (x). From the definition 

[/„(*) = 2" • TT z-cos-^- 

tA\ n+1 

it follows that Ui m {x) = (— l) m ■ i?m (— 4x 2 ). From there, a simple derivation gives 
the generating function of the polynomials Hm\x). 

Proposition 1. The generating function of the sequence Hm\x) is rational: 

V H^(x)t m = 1 — 

^ " w {\-ty-xf 
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Proof. Starting from the classical rational generating function 

f 



(1.3) j2 Un ^)t n 



l-2xt + t 2 ' 

a generating function for the even part is readily obtained: 

V U 2 (x)t^ - ( l - + I ^ - -±l 

^ 2ml 1 ~ 2 I l-2xt + i 2 + l + 2xt + P J (l + i 2 ) 2 -4x 2 t 2 ' 

m>0 v / \ J 

and thus we deduce 

m>0 m>0 ^ ' 

from which the conclusion follows. □ 

In [6] , Proposition f was proved in a different way, and the generating function 
for the case s = 2 was shown to be rational. We generalize these results by proving: 

Theorem 1. For any s > 1, the generating function F s (x,t) = X)m>o Hm\x)t m 
belongs to <Q(x,t). Its denominator has degree at most 2 s in t, and its numerator 
has degree at most 2 s — 1 in t. 

Moreover, we give an algorithm that computes these rational functions explicitly. 
In Section 2, we prove this theorem. In Section 3, we comment on the computational 
aspects. We conclude with a few further observations in Section 4. 



2. Rational Generating Series 

For convenience, we work with the polynomial G$(x) = {-l) m H$ (-x). It 
is monic, with roots the s-th powers of the roots of Gm\x) = {—l) m H^}{—x). 
To prove the theorem, it is clearly sufficient to show that the generating function 
J2 m >o G ™ ( x ) tm is rational. 

2.1. Roots. Letting e s be a primitive s-th root of unity, we have 

g^v)= n (* s -« s )= n (n(*-a/4) ) > 

G^(a)=0 GS'(a)=0 V=° / 

which equals 

ri( n (*-«/^.)] =n^-n G s ) (^)=(-i) mM -nG«(ejx). 

3=° \G, ( ^ , (a)=0 / 3=0 3=0 3=0 

Therefore, the polynomial Gm can be expressed using solely Gm by 

s-1 

(2.1) G&V) = (-l)-(-D • G$(eix), 



3=0 



and in particular it is given by the resultant Gm\x) = Res y (Gm\y), x — y") 
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2.2. Polynomials, Recurrences and Hadamard Products. A consequence 
of Proposition 1 is that the sequence of polynomials Gm (x) satisfies the linear 
recurrence 

(2.2) G < i\ 2 (x) + {2-x)G^ +l {x)+G^{x) = Q, G (x) = 1, Gx(x) = x - 1. 

The coefficients of this recurrence are constant, in the sense that they do not depend 
on the index m. The characteristic polynomial of the sequence Gm (x) is the 
reciprocal i 2 + (2 — x)t + 1 of the denominator of its generating function. 

Conversely, any sequence satisfying a linear recurrence with constant coefficients 
admits a rational generating function and this applies in particular to the sequence 
of polynomials Gm\e J s x) for any j > 0. 

The product of sequences u n and v n that are solutions of linear recurrences with 
constant coefficients satisfies a linear recurrence with constant coefficients again 
(see, e.g., [7, §2.4], [4, Prop. 4.2.5], [2, §2, Ex. 5]). Its generating function, called 
the Hadamard product of those of u n and v n , is therefore rational. From Eq. (2.1), 
it follows that the generating function J2m>o ^rn {x s )t m is also rational (a priori 
belonging to C(x,t j). 

Moreover, the reciprocal of the denominator of a Hadamard product of rational 
series has for roots the pairwise products of those of the individual series. Thus, let- 
ting a±(x), oi2(x ) denote the roots of (1 + t) 2 —xt, the reciprocal of the denominator 
of the generating function of (— l) m ( s_1 ) ■ Gm\x s ) is the characteristic polynomial 

P s (x, t) = (t- a n (x)a i2 (e s x) ■ ■ ■ a ia (e^xfj . 

l<»i,...,i,<2 

Note that, since a\ • a.% = 1, the polynomial P s is self-reciprocal with respect to t. 

We prove that the polynomial P s {x, t) belongs to Q[x s , t] by showing that all the 
(Newton) powcrsums of the roots of P s (x,t) belong to For any I £ N, the 

£-th powcrsum is equal to the product 

(2.3) (a{(x) + a e 2 (x)) (a{(e s x) + a e 2 (e s x)) ■ ■ ■ {a^e^x) + aiie^x)) , 

and thus has the form Tg tS (x) = Qi(x)Qi(s s x) ■ ■ ■ Qi(e a s ~ 1 x) for some polynomial 
Qe(x) € Q[x]. The polynomial Ti tS (x) being left unchanged under replacing x by 
e s x, it belongs to Q[x s ], and so do all the coefficients of P s . This can also be seen 
from the expression of T^ )S as a resultant: Ti, s (x) = ReSj,(?/ s — x s , Qe(y))- 

In conclusion, Gm(x s ) satisfies a recurrence with coefficients that are polynomi- 
als in Q[a; s ], thus the series J2 m >o Gm\x s )t m is rational and belongs to Q(x s ,t), 
and thus F s (x, t) belongs to Q(x, t). The assertion on the degree in t of (the denom- 
inator and numerator of) F s (x, t) follows from the form of P s (x, t). This concludes 
the proof of Theorem 1 . 

3. Algorithm 

The proof of Theorem 1 is actually effective, and therefore it can be used to gen- 
erate, for specific values of s, the corresponding rational function F s , in a systematic 
and unified way. 
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3.1. New Proof for the Case s = 2. Wc illustrate these ideas in the simplest 
case s = 2, for which the computations can be done by hand. 

We compute the generating function of k m {x) :— iJ™'(— £ 2 ) = Gm\x)Gm\—x) 
from which that of the sequence Hm\x) is easily deduced. By Proposition 1, the 
generating functions of the sequences Gm\x) and Gm\—x) have for denominators 

gi{t) := (1 + t) 2 - xt, g 2 (t) := (1 + if + xt, 

which are self-reciprocal. The polynomial whose roots are the products of the roots 
of gi and 172 is obtained by a resultant computation: 

Res u ( gi (u),u 2 g 2 (t/u)) = {t-lf+x 2 t{l+tf = l+{x 2 -4)t+(2x 2 +6)t 2 +(x 2 -4:)t 3 +t 4 

It follows that the sequence k m (x) satisfies the fourth order recurrence 

k m+4 + (x 2 - 4)fc m+3 + (2x 2 + 6)k m+2 + (x 2 - 4)fc TO+ i + k m = 0. 

The initial conditions can be determined separately 

fc = 1, ki = 1 - x 2 , k 2 = 1 - 7x 2 +x 4 , k 3 = 1 - 26x 2 + 13x 4 - a; 6 

and used to compute the numerator of the generating function of the sequence of 
(2) 

polynomials H m (x). In conclusion, we have just proven the identity 

<iU) J^'-' ^-^X+ir 

This provides a simple alternative proof of [6, Prop. 4.1]. 



3.2. General Algorithm. The computation of Hadamard products extends to 
linear recurrences that can even have polynomial coefficients. It is implemented in 
the Maple package gfun [3]. This makes it possible to write a first implementation 
that produces the generating function of the Hm\x) for small s. 

The computation can be made more efficient by using recent algorithms specific 
to linear recurrences with constant coefficients [1] (or, cquivalcntly, to rational func- 
tions). The central idea is to avoid the use of a general-purpose algorithm for the 
computation of bivariate resultants and recover the denominator from the (New- 
ton) powersums of its roots, themselves obtained from the powersums of a\ and ca- 
using very simple univariate resultants, computations in algebraic extensions can 
be avoided. This is summarized in the following algorithm that follows the steps of 
the constructive proof of Thm. 1, using efficient algorithmic tools. 

Input: an integer s > 2. 

Output: the rational function F s (x,t). 

(1) Compute the powersums Qt{x) := ai(x) e + a- 2 {xY e Q[x] for £ = 0, . . . , 2 s ; 

(2) Compute the powersums X^ s from Eq. (2.3) by Zi sfx 1 / 8 ) = RcSj / (y s — 
x,Q e (y)) for £ = 0,..., 2 s ; 

(3) Recover the polynomial P s {x 1 / S ,t) from its powersums T^ s [x x / S , t); 

(4) Deduce the denominator D s (x, t) of F s (x, t) using P s (x 1 ^ s , t) = D s (—x, (— l) s 

(5) Compute G 1 ^ {x) for m = 0, . . . , 2 s - 1 using the 2nd order recurrence (2.2); 

(6) Compute G$(x) = Res y {G^ ] (y), x - y s ) for m = 0, . . . , 2 s - 1; 

(7) Compute the numerator N s (x, t) := A>(M) x Em=o G^{-x){-t) m modt r 

(8) Return the rational function F s (x,t) = N s (x,t)/ D s (x,t). 



CHEBYSHEV-LIKE POLYNOMIALS 



5 



Steps (1) and (3) can be performed efficiently using the algorithms in [1]. In Step (4) 
we use the fact that P s (x,t) is the denominator of the generating series 

. G ^) {x s )t m = j2 H$(-x s ){(-iyt) m = f s (-x°, (-iyt). 

m>0 m>0 

We give the complete (and remarkably short) Maple code in an Appendix. 

3.3. Special Cases. For s = 2, the computation recovers (3.1). For s = 3, it takes 
less than one hundredth of a second of computation on a personal laptop to prove 
the following result 

V „ (3) , v, ffl = (l-f)((t-l) 6 -^ 2 (t + 3)(3t+l)) 

m W ^ _ ^4 + 14t 3 + 34t 2 + Ut + ftU _ ]\2 + ( t _ J)8 ' 

m>0 ' 

that was conjectured by Stolarsky in [5]. 

Similarly, less than two hundredth of a second of computation is enough to 
discover and prove the following new result. 

Proposition 2. The generating function Em>o^""'( ;c )' 1 ™ * s equal to 

(t - 1) (x 2 t 4 A(t) - 2xt 2 B(t) (t - If + {t- l) 14 ) 
x 3 t 5 (t + l) 2 (t - l) 4 + x 2 t 3 C(t) + xt(t - l) 8 D(t) -(t- l) 16 ' 

where 

A(t) = 9t 6 - 46i 5 - 89i 4 - 260t 3 - 89i 2 - 46f + 9, 
B(t) = lit 4 + 128t 3 + 266t 2 + 128* + 11, 

C(t) = 2t 10 - 13t 9 + 226t 8 - 30(M 7 - 676t 6 - 2574t 5 - 676t 4 - 300t 3 + 226t 2 - 13t + 2, 
D(t) =t 6 + 60t 5 + 519t 4 + 1016t 3 + 519t 2 + 60i + 1. 

The code given in the appendix makes it possible to compute all the rational 
functions F s (x,t) for 1 < s < 7 in less than a minute. 



4. Further Remarks 

( s) 

While computing with the polynomials H, n wc experimentally discovered the 
following amusing facts. 



Fact 1. For all s > 1, we have: 

4 s \x) = l, 

H[ s \x)=x+1, 

H^ s \x) = x 2 + L 2s x + l, 

where L„ denotes the n-th element of the Lucas sequence defined by Lq = 2, L\ = 1 
and L n = L n -\ + L„_2 for all n > 2. 

This is easy to prove. 
Fact 2. Hm\x) has only non-negative integer coefficients, for all m > and s > 1. 



(i 
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Non- negativity is clear from (1.2), while the integrality follows from the generat- 
ing function when s = 1 and from the resultant representation of Hm for higher s. 

Fact 2 suggests that the coefficients of the polynomials Hm\x) could admit a 
nice combinatorial interpretation. 

For instance, the coefficient of x 1 in H$ (x) is equal to the trace of the matrix 
M 2s , where M = (a 2J )™ =1 is the m x m matrix with o»j = 1 for i + j < m + 1 
and Ojj = otherwise. 

Fact 3. The rational function F s (x,t) = X) m >o (x)t m writes N s (x,t)/D s (x 1 t), 
where N s (x,t) is a polynomial in Q[x,t] of degree 2 s — 1 in t and B s -\ — 1 in x, 
D s (x,t) is a polynomial in Q[x,t] of degree 2 s in t and -B s -i in x, and where B n 
is the central binomial coefficient 

Here is the sketch of a proof. We only prove upper bounds, but a finer analysis 
could lead to the exact degrees. 

First, it is enough to show that the degree in x of the polynomial P s (x, t) defined 
in the proof of Theorem 1 is at most sB s -\. Now, instead of considering P Sl we 
study the simpler polynomial 

C s (x,t) = Y\ (t-a ll (x)a i2 (x)---a u (x)y 

l<ii,...,i s <2 

We will prove the bound s£? s _i on the degree of C s in x. Adapting the argument 
to the case of P s is not difficult. 

The starting point is that, when x tends to infinity, a\(x) grows like x, while 
ct2{x) grows like x . Therefore, 

C s (x,t) = f[(t-a s 1 - k (x)al(x)y k) 

k=Q 

grows like L s (x,t) = YV k=0 (t — x s ~ k (x~ 1 ) k ^ k ^ . The degree in x of C s (x,t) is thus 
bounded by the degree in x of the Laurent polynomial L s {x,t), which is equal to 

deg x (i s ) = £ Q (s - 2k) = ( [s/2 ] + x ) • (L./2J + 1) = sB s _ v 

Appendix: Maple Code 

For completeness, we give a self-contained Maple implementation that produces 
the rational form of the generating function F s (x,t). 

recipoly :=proc(p,t) expand(t~degree(p,t)*subs(t=l/t,p)) end: 

newtonsums : =proc (p , t , ord) 
local pol, dpol; 

pol : =p/lcoef f (p ,t) ; 

pol : =recipoly (pol ,t) ; 

dpol : =recipoly (dif f (pol ,t) ,t) ; 
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map (expand , series (dpol/pol , t , ord) ) 
end: 

invnewtonsums : =proc (S , t , ord) series (exp (Int ( (coef f (S,t,0)-S)/t,t)) ,t, ord) end : 
Fs:=proc(s,x,t) 

local H, G, N, Q, T, ell, y, Ps, Ds, Gl, Gs , m, Ns; 
H:=(l-t)/((l-t)~2-x*t) ; 
G:=subs(x=-x,t=-t,H) ; 
N:=2~s; 

Q:=subs(x=y,newtonsums(denom(G) ,t,N+D) ; 

T:=series(add(resultant(y~s-x,coef f (Q,t,ell) ,y)*t~ell,ell=0. . N) ,t,N+l) ; 
Ps:=convert(invnewtonsums(T,t,N+l) ,polynom) ; 
Ds:=subs(t=(-l)~s*t,x=-x,Ps) ; # denominator 
Gl : =subs (x=y, series (G, t ,N) ) ; 

Gs : =series (add (resultant (coef f (Gl ,t ,m) ,x-y~s ,y) *t~m,m=0 . .N) ,t,N) ; 
Ns:=convert (series (Ds*subs(x=-x,t=-t,Gs) ,t,N) ,polynom) ; # numerator 
collect (Ns/Ds,x, factor) 
end: 
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